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INDCCDCXKM 


Progress by investigators at the Georgia Institute of Technology in the 
development of techniques for passive microwave retrieval of water vapor 
and precipitation paramet ers using millimeter- and sub-mil 1 imeter 
wavelength channels is reviewed. Channels of particular interest are in the 
tropospheric tra nsmiss ion windows at 90, 166, 220, 340 and 410 GHz and 
centered around the water vapor lines at 183 and 325 GHz. Collectively, 
these channels have poten tial application in high— resolution mapping / 
from geosynchronous orbit) , remote sensing of cloud and precipitation 
parameters, and retrieval of water vapor profiles. 

During the period from January 1, 1992 through June 30, 1992, the 
Millimeter-wave Imaging Radiometer completed successful integration and 
data flights on the NASA ER-2. Coincident data was collected from several 
other sources, including the EMSP SSIVT-2 satellite, a collocated 
radiosonde, the MIT Millimeter-wave Temperature Sounder (MIS) and the 
na<*ffaing infrared imager. The mechanical and electrical integrity of the 
MIR are quite good, and the sensitivity for the six functioning channels 
(90, 150, 183+/-1,3,7, and 220 GHz), is close to specification. Integration 
of the 325+/-l,3, and 9 GHz channels should be completed during the fall of 
1992 during preparation of the MIR for TOGA/GQARE. 

In order to resolve q ues tions cor*.*-: riling the abs olute calibration of 
the MIR, reduction of laboratory bistatic scattering measurements of 
microwave blackbody calibration loads has continued. Favorable comparisons 
between the reduced data for one-dimensionally periodic structures and 
predicted results using the coupled-wave method suggest that the coupled 
wave method will be useful for characterizing both one- and two- 
dimensicnally periodic loads. A comprehensive thermal analysis is still 
needed. 

New studies commenced under this grant include the investigation of the 
utility of the third Stokes parameter Ty=<E v I^ 1 *> for passive remote sensing 
of the surface state. Here, we have shewn using laboratory wave tank 
measurements that a disoemable signal (up to +/— 10 K) related to water 
wave direction can be obtained using a polarization correlating radiometer. 
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SOMARY OF ACTIVITIES 


Activities within the period frtxn January 1, 1992 through June 30, 1992 
by Georgia Tech researchers in millimeter and submillimeter wavelength 
tropospheric remote sensing have been centered around the integration and 
initial data flights of the MIR on board the NASA ER-2. Georgia Tech 
contributions during this period include completion of the MIR flight 
software and implementation of a "quick-view” graphics program for ground- 
based calibration and analysis of the MIR imagery. 

In the current configuration, the MIR has channels at 90, 150, 183+/- 
1,3,7, and 220 GHz. Provisions for three ad d it ional channels at 325+/-l,3 
and 9 GHz have been made, and a 325-GHz receiver is currently being built 
by the ZAX Millimeter Wave Corporation for use in the MIR. The combination 
of the mil l imeter wave and submillimeter wave channels aboard a single 
well -cal ibrated instrument will provide the necessary aircraft radiometric 
data for radiative transfer and cloud and water vapor retrieval studies. A 
paper by the PI discuss ing the potential benefits of passive millimeter- 
and s ubmi 1 1 wave observations for cloud, water vapor and 
precipitation measurements has recently been acc e pted for publication 
([Gasiewski, 1992], see Appendix A) . 

The mtr instrument is a joint project between NASA/GSFC and Georgia 
Tech. Other Georgia Tech contributions to the MIR and its related 
scientific uses have included basic system design studies, performance 
analyses, and circuit and radiometric load design. 

1. MIR Flight and Quick View Software 

At this time, all essential flight software for unattended operation of 
the MIR aboard the NASA ER-2 has been completed. The robustness of the 
flight software has been demonstrated by three data flights during May, 
1992, totaling 13 flight hours. During the first flight (the integration 
test flight) of the MIR, seme erratic mirror movement occurred using the 
earliest version of the fli^it software, rendering only approximately 90% 
of the data useful. However, only minor debugging was necessa r y to correct 
this, and no abnormal operation was noted during two successive data 
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flights. 

A summary of the May 1992 flights is shewn in Table 1. Applicable data 
from other instruments, including the Massachusetts Institute of 
Technology's Millimeter-wave Temperature Sounder (MIS), the NASA Airborne 
Ocean Color Imager (ACd) , the EMSP SStyT-2 radiometric imager, and 
collocated radiosondes is available and currently being compiled. 

To expedite post- flight data analysis, a "quick-view" graphics program 
for data display and analysis was written. The program reads flight data 
tapes and displays the brightness imagery in six 61x288 pixel maps, 
comprising approximately 28 minutes of data. Brightness maps appear with a 
nearly 1:1 aspect ratio using a 16-oolor scale. A variety of data editing 
features are available. 

Calibration is accomplished using hot and cold load views (once per 
scan) along with appropriately weighted temperature values from the eight 
FID blackbody load sensors. A causal infinite impulse response (IIR) filter 
is currently u sed to estimate instrument gain and offset values. The 
software is being upgraded to include an op timal non-causal Wiener filter 
for gain and offset estimation. Either absolute brightness temperatures or 
brightness temperature structure can be displayed using the program. The 
structure display is particularly useful for analyzing small perturbations 
in bright ness caused by clouds or water vapor f eat ures. 

Typical brightness structure imagery far the 89, 150 and 220 GHz 
channels from sortie 92-090 (May 16, 1992 , 0042-0059 UTC) is shewn in 
Figure 1. This imagery was observed over the Pacific Ocean west of Los 
Angeles, CA, and is displayed at approximately the Nyquist sampling rate 
for the 3.5° antenna beamwidth of the MIR. Ihe imagery reveals striated 
cloud structures as gmai l as a few kilometers in size (A) , and numerous 
small lcw-level clouds (B-E) . A small island (San Nicholas) is seen in all 
channels. Brightness perturbations from the nominal baseline temperatures 
are as shown in Table 2. Much of the low cloud structure seen in the 89 and 
150'GHz channels is screened by water vapor, as evidenced by the lack of 
structure in the 220-Giz channel. Beyond the clouds (near the right of the 
images) , a marked decrease in 220-GHz brightness suggests a corresponding 
reduction in water vapor. 
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2. W-band load Sc attering Measurements 

During the fall of 1991, bistatic scattering measurements at W-band 
were mark* for several microwave absorbing surfaces, i ncludin g flat, wedge- 
type, and pyramid-type iron-epaxy blackbody absorbers and flat carbon-foam 
"Echosorb" samples. Hie measurements were made using a Hewlett-Packard 8510 
network analyzer interfaced to a focussed-lens scattering range at the 
Georgia Tech Research Institute. Hie dynamic range of this facility is 
approximately 55 dB, and the accuracy of each measurement is estimated to 
be 1 to 2 dB. Swept- frequency measurements over the range 75-100 GHz have 
revealed and Bragg reflection characteristics for the wedge-type 
sairples [Gasiewski and Jackson, 1992 (Appendix B) ] . 

Hie pr imar y purpose of the bistatic scattering data is to verify the 
electromagnetic scattering calculations obtained using the coupled wave 
method [Gasiewski et al, 1991; Moharam and Gaylord, 1982]. In order to use 
the laboratory data, deconvolution of the spatial response of the antennas 
used in the scattering range was required. Hiis was performed using the 
CLEAN algorithm [Thompson, et al, 1991], which is particularly applicable 
for deconvolution of blurred point-source imagery. In the periodic absorber 
case , the Bragg harmonics exhibit scattering peaks as a function of angle 
whirh can be tre ated in an analogous fashion to point sources. A comparison 
of the deconvolved laboratory data and coupled wave calculations is 
illustrated in Figure 2. Here, the calculated and measured reflected powers 
in both the gpomiar- and first-order Bragg harmonics are seen to agree to 
within 5-dB, which is approximately the accuracy to be expected from this 
measurement. The favorable agreement validates the coupled wave method and 
suggests that its extension to two dimensions will permit detailed analysis 
of pyramidal blackbody loads. 

Unfortunately, scattering from cur single pyramid load sanple 1 is not 
measurable using the current range configuration. This is not to say that 
the pyramid load does not scatter, but rather that the scattered energy is 
distributed over a large range of solid angle. Thus, a smaller flux density 
is typically presen t in any given scattering direction, neces s itating a 


1 The pyramid load sanple has been provided by the ZAX Millimeter-wave 
Corporation, San Dimas, CA. 
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more sensitive receiver. Future measurements will use a radiometric 
receiver in place of the HP 8510 W-band dcwnocnverter. This should improve 
the sensitivity of the measurement by 35-40 dB. 

3. 325-GHz Receiver for use on the NA5A/GSCF MIR 

The implementation of a 325/340 GHz receiver on the MER has been a high 
Georgia Tech priority. A double-balanced Schottky diode mixer employing a 
subharmcnically-punped local oscillator Gunn diode source is currently 
being fabricated by the ZAX Millimeter-Wave corporation of San Dimas, CA. 
The double-sideband noise figure quoted far the receiver is 12 dB (4300 K 
noise temperature) , which is acceptable for aso e r t a i nin g science results at 
S3MMW freq ue ncies. The receiver should be available for testing and 
integration into the MIR near the end of the summer of 1992, and operative 
for ER-2 data flights during T0GA/CQARE, January-February, 1993. We are 
currently studying the design of a plasma noise source for testing the 325- 
GHz receiver. 

4. Polarized Microwave Radicmetrv 

Additional experiments that we have performed under funding from this 
grant include measurements of microwave emission from a striated water 
s u rta x using a 90-GIz tri-polarizaticn radio me te r . The radiometer, 
developed under a previous NASA grant (NAG 8-829), permits precise 
measurement of the first three Stokes' parameters. This is accomplis h ed 
using a conventional dual-polarization radiometer augmented with a third 
cross-correlation channel. Accurate calibration is performed using an 
innovative polarized calibration load [Gasiewski and Kunkee, 1992a]. Using 
the three channels, the polarization basis can be arbitrarily rotated to 
eliminate the polarization skew that accompanies a mechanically-scanned 
imager. This technique is termed electronic polarization basis rotation 
(EFBR). 

While the utility of vertically and horizontally polarized brightness 
temperatures for both surface and atmospheric remote sensing has long been 
recognized, the potential use of the third Stokes parameter l\j‘=<E v Ef 1 *> in 
Earth remote sensing has largely been overlooked. However, our recent 
laboratory measurements vising a fresh-water wave tank illustrate the 
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dependence of T^j on the angle of the water wave (Figure 3) . The T^j 
experimental data for the fas** of a 65° observation angle is well 
corroborated using the Kirchoff model for surface reflectivity, although 
not all incident angles and polarizations agree quite as well. 

The relatively large anplitude of the Tfy signal (+/-10 K) suggests that 
remote sensing of ocean wave direction might be facilitated using a 
spacebome tri-polarization (as opposed to dual polarization) radiometer. 
In addi tion, the measurement of T^j (along with T v and T^) facilitates EPBR, 
thus eliminating the need for the space-qualified bearings and slip rings 
that are used on the EMSP SSty/I imaging radiom e ters. 
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GGNCHSICNS AND FIAN5 PCR FUTURE WCRK 


To date, Georgia Tech has provided assistance to NASA/GSFC in the 
mechanical and electrical design of the Mik, and has been pr imar ily 
responsible for MIR flight and analysis software. At this ti m e, cur pr i m a r y 
interests are turning toward the integration of the 325-GHz receiver into 
the MIR, airborne data collection using the MIR, and the development of EHF 
and SMMW meteorological retrieval algorithms . The radiometric calibration 
load reflectivity analysis will be continued, although this study will be 
of tertiary importance relative to the MM and SMMW data collection and 
analysis. The highest priority is assigned to procuring and integrating the 
325-GHz receiver into the MIR immediately after planned aircraft flights 
during the summer of 1992. At this time, it appears hi^ily likely that the 
325-GHz channels will be operational for TOGVOQARE in early 1993. 

Concurrent with the MIR operation, retrievals of water vapor and cloud 
water content fr om EHF and SMMW observations will be attempted. The most 
premising method seems to be the non— li near sta t isti c al iterative method 
used by Kuo [1988]. An extension of this technique to include cloud water 
content, base altitude, cloud thickness and surface eraissivity as variable 
parameters is currently being studied. To provide sufficient ocrputational 
capabilities (disk storage capacity and processing speed) to perform the 
mtr data p rocess ing and retrieval studies, we are planning to purchase a 
Hewlett-Packard model 710 workstation (or equivalent) under an increment to 
this grant. 

To accur a tely calibrate the MIR, the total calibration load 
reflectivity needs to be less than 1%, and known to better than 0.1%. 
Manufacturer's specifications typically provide only the specular cernponent 
of the reflectivity, which is thought to be substantially less than the 
total reflectivity. In order to refine the Mik calibration, we plan to 
extend the study of the electromagnetic characteristics of wedge-type 
blackbody loads to the more desirable pyramidal loads. CXir approach is to 
develop numerical models for one— and two— dimensionally period i c lossy 
gratings using the coupled wave method. 

We currently have software based on the coupled wave method to predict 
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the reflectivity of cne-dimensicrally periodic loads of arbitrary 
dielectric profile, Comparisons of the ocmputed Bragg-lobe reflectivities 
with deconvolved measurements obtained from the focussed-lens scattering 
range show good agreement. Using the ccupled-wave software, optimization of 
the dielectric profile far minimum reflectivity over wide microwave bands 
is envisioned. This has been successfully s h ed for wedge - type 
loads; extension of the coupled wave model to two-dimensionally periodic 
surfaces appears feasible, and will be attempted. 

oi or^rrmag net.i c analysis of calibration l oads has al so benefited 
other involved in micrcwave radicmetry. Par example, cur numerical 
program was recently (September, 1991) used to guide engineers at NASA/MSPC 
in the selection of load geometry for new wideband (10-85 GHz) calibration 
loads for the Advanced Micrcwave Precipitation Radiometer (AMER) . Future 
plans include a study of the steady-state thermal profile in both wedge- 
and pyramid-type loads for the purpose of predicting the emitted thermal 
radiation. 

Tto further investigate the potential for tri-polarization radicmetry in 
the r'° n r ri ~ e ‘ sensing of surface and at nfn Aftr ic par ameter s , three additional 
projects are planned; (1) Ground-based sky observations of thunderstorm 
anvils are to be made during the months of July and August, 1992 using the 
tri-polarization radiometer. Of particular interest will be correlations 
between polar imet ric signatures in the micrcwave data and thunderstorm 
electrification. A field mill probe is to be supplied by NASA/MSPC for 
correlation purposes. (2) To reduce the complexity of the hardware 
accnri a+oH with calibrating a tri-polarization radiometer, a digital cross- 
correlator is being designed. The cross correlator will be able to be used 
on NASA DC-8 ai r craft flights of the 90-GHZ tri-polarization radiometer, 
and will be a precursor to the one proposed to be used on the NASA/MSPC 
Advanced Micrcwave Precipitation Radiometer (AMER) [Gasiewski and Kunkee, 
1992b]. The bandwidth of the cross-correlator will be approximately 150 
MHz, makin g the device useful for polarization correlation radicmetry 
radiometric interferometry. Completion of the cross-correlator for use 
during TOGA/GQARE is planned. (3) To assess the potential of tri- 
polarization radicmetry for remote sensing of surfaces and vegetation, the 
surface measurements will be repeated using a striated surface of wet and 
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dry sand in place of the water waives. 

Recommended Improvements to the MIR 

In view of the performance of the MER during aircraft flights and 
ground-based calibration, several minor modifications are recommended to 
improve stability and calibration accuracy. These include: 

1. Transverse thermal gradients of up to 1 K were observed in the hot 
calibration load in flight. These can be reduced by t herm a l ly insulating 
the load f rom its metallic shroud, insulating the shroud using a foam 
blanket, insulating the rear of the load, and finally, standing the load 
off from its mounting pedestals losing t h e rm a l ly insulating spacers. 
Acceptable (and attainable) gradients are appro x i m a t ely 0.1 K. 

2. To minimiz e calibration uncertainty due to load reflectivity and 
spillover, it is recommended that both loads be moved closer to the 
scanning mir ror (approximately 2 cm is available) , the scan cavity be lined 
with moderately absorptive material (e.g. "Echosarfo") and its approximate 
temperature monitored, and the antenna beams (behind the lens array) be 
isolated from each other using thin (2-ram) wadis of absorbing material. 

3. For precise radiative transfer model interocrparisons, the RF 
passband response of each channel (part j rail ariy the water vapor channels 
near 183 and 325 GHz) should be measured. This can be acocrplished using a 
plasma discharge noise source and a spectrum analyzer. 

4. The digital -to-analog offset circuitry should be rewired to take 
full advantage of the eight D/A channels available. This would req u ire , for 
example, coupling the 183+/— 1 and 183+/— 3 channels, with all other channels 
having independent offset lines. A minor flight software modification would 
be required to aocxximodate the wiring change. 

Other slightly more difficult modifications and/or tests include: 

5. Pre ci se in-situ sensitivity patterns for all receivers should be 
measured, i ncludin g both near and far— zone sensitivities. This will allow 
the determination of the amount and significance of beam spillover at the 
calibration loads. If needed, new load shrouds might be machined to subtend 
a slightly larger solid angle. 

6. A grooved TEX insulator should be designed and fabricated to cover 
the calibration loads. Using the coupled-wave software, the transparency of 
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the insulator could be optimized for the suite of MIR frequencies. It is 
expected that such a cover w ill exhibit less scattering and absorption than 
the s t y rofoam ocrvers currently used, yet provide adequate t h e rm a l 
insulation. Provisions for purging the space behind the load cover with dry 
nitrogen before flights should be installed. 

7. During flights, sane instability in the temperatures of the 
calibration l oads and miver temperatures was noted. Par this reason alone, 
the temperature regulators should be redesigned the nonlinear dynamic 
thermal responses of the various thermo-electrical circuits considered. 
However, the current regulators are also quite inefficient in that they are 
resistive (lossy) in design. Hence, they heat the MIR frame im m ediately 
surrounding them considerably. Modified regulators based on lossless 
switching designs should be used here. Provided proper grounding and 
shielding practices are followed, interference caused by switching noise 
from these devices should not occur. 
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Table 1. 


MIR integration and data flights: Spring 1992. 

Sortie # Date Time Instruments 

(17TC) 

fwe) 

92-087 5/11/92 1900-2200 MIR, MIS (u) ,AOCJ 

92-089 5/14/92 2315-0515 MIR, MIS (d) ,AOCX 

92-090 5/15/92 2320-0500 MIR,MES(u) , A0CX 


Norte: "u" or "d M indicate up-looking or down- 
looking, respectively. 
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Table 2 


Brightness perturbations (in K) far various features in Fig. 1. 


89 150 220 

location GHz GHz (3iz 


A 

20 

8 

0.5 


B 

25 

14.5 

1 


C 

22 

14.4 

1 


D 

18 

11 

0.7 


E 

23 

13 

0.7 


I (island) 

67 

25 

5 


Baseline (nadir) 

208.5 

255 

279 
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~ 190 km 


<r 




A 



89 GHz - range=25 K 


* 


~ 40 km 
at 

surface 


* 



150 GHz - range =25 K 



220 GHz - range =7 K 


Figure 1: Brightness structure maps observed from 0042 to 0059 
OTC, May 16, 1992 over the Pacific Ocean during sortie 92-090. 
The letters designate brightness features listed in Table 2. 
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Reflected Power (dB) 


Wedge Absorber (J) 
TE 60° Incidence 



Figure 2. Comparison of the measured reflected power at 89 (Hz 
from a one -dimen sionally periodic blackbody absorber with wedge 
cross-section. The incident mode is TE with a 60° incidence angle 
with the grating vector in the plane of incidence. Ihe load 
sample is 6.5 mn deep, 3.2 mm in pitch, and the relative 
permittivity and permeability are approximately 9-j0.4 and l-j0.5 
(respectively) . 
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3 rd Stokes Parameter Ty vs Wave Angle 

1 — cm fresh water standing waves: A=10 cm 
f = 92 GHz, 65° incident angle 



Fiaure 3: Ccrpariscn of laboratory measuranents of the third 
LokS peu^S^cteerved over a striated fresh watersurf ace 
with calSlations usi ng the Kirchoff approoomation. The abscissa 
is the azimuthal angle of the water wave direction. The angle of 
observation is 65° from the nadir. 
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Abstract 


Potential uses of specific Extremely High Frequency (EHF) and Sub- Millimeter- 
Wave (SMMW) channels at 90, 166, 183, 220, 325, 340, and 410 GHz for passive 
spaceborne remote sensing of the troposphere and lower stratosphere are investigated 
using an iterative numerical radiative transfer model. Collectively, these channels 
offer potential for high spatial resolution imaging using diffraction-limited apertures 
of practical size, along with the ability to profile water vapor, map precipitation 
beneath optically opaque cloud cover, and to measure nonprecipitating cloud (e.g., 
cirrus) parameters. The numerical study, along with the results of previous passive 
microwave experiments, suggests that a widely-spaced set of EHF and SMMW chan- 
nels can yield observable degrees of freedom related to clouds and precipitation not 
available by exclusively using the more thoroughly studied microwave channels below 
183 GHz. A new passive airborne imaging instrument for tropospheric meteorolog- 
ical sensing at 90, 150, 183±1,3,7, 220, and 325±1,3,9 GHz, the Millimeter-wave 
Imaging Radiometer (MIR), is described. 
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1 Introduction 


Spaceborne passive microwave measurements of tropospheric meteorological 
parameters are impeded by at least two problems: (1) spatial undersampling of the 
brightness field due to the inherently broad antenna patterns of practical satellite- 
based microwave imagers, and (2) the adverse effects of heavy, non-precipitating 
clouds on both temperature and water vapor retrievals. On occasion, it is also de- 
sirable to have enhanced sensitivity to thin clouds (e.g., cirrus) for the purposes of 
detection and water content estimation. This need contrasts with the primary advan- 
tage of microwave observations in relation to infrared (IR) and visible observations, 
namely, an enhanced ability to probe through most cloud cover. 

These problems can be partially alleviated by observing at specific Extremely 
High Frequency (EHF: 30-300 GHz) and Sub- Millimeter Wave (SMMW: 300+ GHz) 
passive microwave frequencies. Of particular interest are the EHF and SMMW chan- 
nels at 220, 325, 340, and 410 GHz, used in conjunction with the more thoroughly 
investigated (or, “conventional”) microwave channels at 6, 10, 18, 37, 50-60, 90, 118, 
166, and 183 GHz. Collectively, these channels sample nearly all significant clear-air 
tropospheric spectral features (i.e., all O 2 and H 2 O absorption lines and transmission 
windows, see Fig. 1). In particular they sample the relatively featureless hydrometeor 
absorption and scattering spectrum over nearly two decades in frequency, spanning 
the transition from the Rayleigh to the Mie regions for a wide range of precipitation 
rates (Figs. 2a and b). Interest in EHF and SMMW channels stems from both (1) the 
increased spatial resolution available using diffraction-limited apertures of fixed size, 
and (2) hypothesized observable degrees of freedom related to water vapor, clouds, 
and precipitation parameters using these channels both alone and coincident with 
other lower-frequency microwave channels. 

This paper discusses the results of a numerical investigation into the tropo- 
spheric meteorological observables that can be expected using 220, 325, 340, and 410 
GHz satellite-based brightness measurements, both alone and in coincidence with the 
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conventional microwave channels. For purposes of passive sensing, tropospheric ab- 
sorption is negligible below 6 GHz and generally excessive above 450 GHz. Thus, 
channels outside this range are not considered. The following potential capabilities 
of the collective set of channels are addressed: 

1. Imaging with enhanced spatial resolution, 

2. Water vapor profiling using 325 GHz spectra, 

3. Cloud detection and water content estimation, 

4. Precipitation mapping, 

5. Observation of additional meteorological degrees of freedom using a wideband 
channel set. 

A new cross-track scanning imaging radiometer for use on the NASA ER-2 and 
DC-8 aircraft, the Millimeter-wave Imaging Radiometer (MIR), is described. Data 
from the MIR will be used to verify EHF and SMMW radiative transfer models and 
demonstrate tropospheric water vapor and cloud parameter retrieval and precipitation 
mapping using EHF and SMMW channels. 


2 Spatial Resolution Enhancements 


Because of the the non-linear relationship between hydrometeor parameters 
(e.g., rain rate) and brightness temperatures for frequencies above ~ 10 GHz, spatially 
undersampled imagery can result in ambiguous precipitation parameter retrieval. For 
example, at 37-GHz, a 50-mm/hr oceanic rain cell occupying 10 percent of a beam 
footprint produces approximately the same observed brightness as a 1 -mm/hr driz- 
zle occupying the entire footprint (as can be inferred from [3]), even though the 
footprint-averaged rain rates differ by a factor of five. Even in the case of passive mi- 
crowave temperature sounding (which uses a linear retrieval operator), the presence 
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of small undetected precipitation cells can cause substantial brightness perturbations 
that are indistinguishable from mesoscale temperature fluctuations. If the fractional 
coverage of precipitation within the footprint is unknown, unacceptable temperature 
profile retrieval errors can result. Similarly, the presence of undetected clouds is ex- 
pected to cause unacceptable water vapor profile retrieval errors [4, 5]. Since clouds 
and rain cells are typically no smaller in horizontal extent than a few kilometers, spa- 
tial resolution comparable to this is ultimately desirable for clear-air screening and 
unambiguous cloud and precipitation measurements. 

The sub-satellite 3-dB footprint size of a diffraction-limited microwave sensor 
with aperture diameter D (in m) is: 

R = 0.38-^j (1) 

where R is the footprint diameter (in km), h is the satellite altitude (in km), / is 
the frequency (in GHz), and a quadratic aperture field taper is assumed. Footprint 
sizes of satellite- based microwave imagers such as the NOAA Microwave Sounding 
Unit (MSU, 110 km at ~ 50 GHz), the Advanced Microwave Sounding Unit (AMSU- 
B, 15 km at 183 GHz), and the Defense Meteorological Satellite Program (DMSP) 
Special Sensor microwave Imager (SSM/I, 12.5 km at 85.5 GHz) are larger than the 
sizes of many clouds and convective precipitation cells. However, from a 1000-km 
altitude low-Earth orbit, acceptable sub-satellite spatial resolutions of R = 3.5 and 
2.3 km at / = 220 and 325 GHz (respectively) could be obtained with a 0.5-m 
aperture. 

High spatial resolution will be particularly difficult to achieve from geosyn- 
chronous orbit using conventional microwave channels [6]. As an example, the size of 
the largest solid aperture-type antenna that can be launched using either the space 
shuttle or a Titan-IV vehicle is ~ 4.4 m. The smallest footprint size available using 
a this aperture at 90 GHz is R = 34.5 km. However, at 325 and 410 GHz these sizes 
are R = 9.5 and 7.6 km (respectively), and are satisfactory for many hydrological and 
meteorological requirements [7]. 
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3 Water Vapor Profiling Using 325 GHz Spectra 


The microwave absorption spectrum of the asymmetric top 1 H 2 16 0 is due to 
rotational transitions induced by the interaction of external fields with the molecule’s 
permanent electric dipole moment. Water vapor resonances at microwave frequencies, 
in particular 183.310 and 325.153 GHz, are adequately modelled by the Van Vleck- 
Weisskopf line shape function [8]. Additional microwave absorption contributions by 
water vapor have been hypothesized to arise from the wings of strong H 2 O rotational 
lines located in the far-infrared. In the microwave region, these contributions form a 
“continuum” that varies slowly with frequency. To date, attempts to model the H 2 O 
continuum absorption using both statistical and impact theories [9] have had limited 
success. However, extensive empirical characterization of both H 2 O continuum and 
resonant absorption has been done by Liebe [10], whose reduced-line water vapor 
absorption model [11, 12] is adequate for frequencies below 1000 GHz: 
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where v is in GHz, and 0 = 300 /T, with T in K. The term in braces is the H 2 0 
contribution (in ppm) to the imaginary portion of the complex air refractivity. The 
parameters 61 *, f» 2t , 63 ,, and center frequencies r /,• for the seven most significant H 2 O 
lines for tropospheric and lower stratospheric passive sensing, along with continuum 
parameters bj and b t for the Liebe model are listed in Table 1 . The parameters of the 
continuum term depend on the particular set of lines that are included in the reduced 
summation. The Liebe water vapor absorption model is consistent with clear-air 
absorption measurements at 337 GHz [13, 14]. 


From the Liebe model, the 183 and 325 GHz lines are seen to have nearly the 
same strength and broadening characteristics. Thus, it is expected that satellite- or 
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aircraft-based nadiral observations near the 325 GHz line can be used to profile water 
vapor in clear-air with essentially the same accuracy as the 183-GHz line. (Profiling in 
the presence of clouds is discussed in Section 6.) This can be seen by considering the 
similarity of computed sets of nadiral clear-air weighting functions for channels near 
183 and 325 GHz. The weighting functions display the response of the * lh channel to 
the kinetic temperature T(h) at various atmospheric levels: 

OO 

T b , = J T(h') Wi(h') dh’ + T s W St (f ) + T cb W CB i(f) (3) 

0 

where T B i is the brightness temperature for the specific passband of the i th channel, 
and the weights and Wcm are the contributions to the observed brightness from 
the surface temperature Ts and cosmic background temperature Tcb — 2.7 K. 

Weighting function computations were made using the Liebe water vapor 
model and the Rosenkranz [15] oxygen absorption model, along with an exponen- 
tial water vapor profile and the U.S. 1976 mean atmosphere temperature profile [16]. 

A 2-km water vapor scale height and land background were assumed. Contribu- 
tions from six SMMW 0 2 lines above 368 GHz (as listed by Liebe [11]) are included. 
Nonresonant N 2 absorption [17] is also included, however, it contributes only a minor 
amount to the total absorption. Contributions from other constituents are practically 
negligible for broadband channels (i.e., bandwidths greater than ~ 10 MHz) below 
410 GHz. Indeed, H 2 0 and 0 2 are the dominant gaseous absorbers in the troposphere 
and lower stratosphere at SMMW and lower frequencies. 

The computations show that temperature weighting functions for double-sideband 
channels near 183rfcl, 3, 7 GHz (Figs. 3a and b) are essentially identical to ones near 
325±1,3,9 GHz, and a channel near the 220 GHz transmission window is essentially 
identical to the 166 GHz channel. Thus, these two channel sets are expected to be 
equally sensitive to water vapor for purposes of clear-air profiling. This equivalence is 
further illustrated by the difference between computed nadiral brightness temperature 
spectra near 183 and 325 GHz for various atmospheric scenarios (Fig. 4). For two 
disparate relative humidity profiles (20% and 80% surface relative humidity, SRH), 
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the clear-air 183- and 325- GHz spectra are offset by a constant spectrum, suggesting 
that the difference spectra are nearly independent of the humidity profile. Hence, it 
is expected that clear-air water vapor profiling can be performed, in principle, equally 
well using channels near either the 183- or 325-GHz line. Similarly, it is expected that 
low- to mid-altitude (~500-800 mbar) water vapor can be probed equally well using 
channels at either the 166- or 220-GHz transmission windows. 

Simulations of the AMSU-B instrument’s capabilities suggest that relative hu- 
midity retrieval errors as low as 5-10% in the 300-1000 mbar pressure region can be 
achieved using single-spot 90, 166, and 183-GHz data [18, 19]. The absence of a low- 
altitude sensing SMMW channel (comparable to 90 GHz) is a slight drawback. This 
will somewhat reduce the 220/325-GHz retrieval accuracy in the the 800-1000 mbar 
pressure region, particularly for summer tropical and mid-latitude scenarios. If de- 
sired, a 90-GHz channel can be included along with the 220/325-GHz channels to 
provide low-altitude sensitivity, albeit at degraded spatial resolution. 


4 Cloud Detection and Water Content Estima- 
tion 


Most non-precipitating clouds consist of hydrometeors of mean volume-equivalent 
radius ~ 50 — 100 /zm and of densities within the range ~ 0.01 to ~ 0.1g/m 3 . Particle 
phases are generally either pure liquid or layered liquid and ice (as in stratus and 
cumulus congestus) or pure ice (as in cirrus). At IR and optical wavelengths, most 
of the constituent hydrometeors have electrical sizes ka » 1, where k = ^ , A is 
the wavelength, and a is the volume-equivalent radius. The associated cloud optical 
depths at IR wavelengths are typically less than one hundred meters [20]. Thus, di- 
rect sub-ceiling probing of non-precipitating clouds using IR channels is not possible. 

In addition, these particles behave as Mie scatterers. Hence, passive IR cloud obser- 
vations exhibit little spectral variation that can be easily linked to parameters of the 
cloud particle distribution. 
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However, for frequencies up to ~ 400 GHz, the electrical sizes of non-precipitating 
cloud particles are typically less than or near unity. Thus, these particles behave 
roughly as Rayleigh scatterers and absorbers, and hence exhibit a monotonically in- 
creasing frequency dependence in both scattering and absorption coefficients. The 
dependence is illustrated in Figs. 5a and b, as based on full Mie calculations of the 
scattering and absorption cofficients for clouds of exponential hydrometeor size dis- 
tribution [21]. For a cloud with a 50 — pm mean particle radius, the liquid extinction 
coefficient increases from ~ 8 dB-km -1 -g -1 -m 3 at 90 GHz to ~ 40 dB-km _1 -g -1 -m 3 
at 410 GHz (a factor of ~ 5). From these values, it is expected that non-precipitating 
liquid clouds will be increasingly detectable at 90, 166, 220, 340 and 410 GHz, yet 
will be transparent to a degree suitable for sub-ceiling measurements of cloud pa- 
rameters. Over the same frequency range, the ice extinction increases from ~ 1 to 
~ 27 dB-km -1 -g -1 -m 3 (a factor of almost 30), suggesting that the SMMW chan- 
nels will be particularly sensitive to ice clouds in comparison to the conventional 
microwave channels. 

The above hypotheses have been investigated by computing the nadiral bright- 
ness temperature changes (relative to clear-air) caused by a 2-km thick cloud of vary- 
ing density, base altitude, latitude, and season. Computations use the iterative radia- 
tive transfer model described by Gasiewski [1], applicable for a planar-stratified Mie 
scattering atmosphere consisting of liquid-, ice-, and mixed-phase clouds. Hydrome- 
teor scattering and absorption are modelled by Marshall-Palmer (MP) and Sekhon- 
Srivastava (SS) distributed [22, 23] liquid or ice phase Mie spheres with Henyey- 
Greenstein (HG) phase functions. Two disparate temperature profiles derived from 
a seasonal and latitudinal interpolation to the 1976 supplemental U.S. standard at- 
mosphere listing are used: January at 90° N latitude (polar winter) and July at 30° 

N latitude (subtropical summer). The relative humidity is nominally 40% at the sur- 
face, with a 2-km exponential scale height. However, within cloud layers, the water 
vapor partial pressure is adjusted to equal the saturation vapor pressure p >a ti the 
determination of which is based on the temperature T and the mass fraction // of 
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total cloud water in the liquid phase: 


Piat(T,ft) = + (1 " fi)P.at,i(T) (4) 

Cloud densities range frona 0.003 to 0.1 g/m 3 and the cloud base-altitude is either 
6 or 8 km. The strong 424.763-GHz 0 2 line cools the 410-GHz window brightness 
temperature by ~ 1 K, and hence is included in the absorption calculation. 

The results of the numerical calculations are shown in Figs. 6a-e, where a land 
background with 5% reflectivity has been used. For polar winter conditions (dark 
circles), all cloud water will be in the ice phase. Under these conditions, the highest 
frequencies considered (340 and 410 GHz) provide the greatest sensitivity to cloud 
water, exhibiting negative perturbations of up to 40 K (solid lines) for 0.02 g/cm 2 
ice integrated water density (IWD). By assuming the ice cloud detection threshold 
to be ~ 1 K (the single-spot sensitivity of a typical satellite-based radiometer), the 
minimum detectable ice IWD is ~ 0.001 g/cm 2 . This is approximately five times 
more sensitive than at 90-GHz. 


However, random brightness perturbations of ±3 — 4 K at 340 or 410 GHz can 
be expected due to variations in atmospheric water vapor content. These variations 
can be parameterized by changes in the surface relative humidity (Fig 6d-e, dashed 
lines, dark circles). In addition, clear-air zenith transmissivities at 340 and 410 GHz 
are 75% and 55% (respectively), hence random variations of ±3% in the surface 
emissivity will produce additional brightness perturbations of ±5 - 6 K. In order 
to unambiguously detect clouds in the presence of both water vapor and surface 
variations, the cloud detection threshold must be increased from ~ 1 K to ~ 8 - 
10 K. In this case, the minimum detectable IWD using 340 or 410 GHz rises to 
~ 0.004 g/cm 2 . At 90 GHz, brightness perturbations due to water vapor variations 
are negligible, but those caused by variations of up to ±3% in surface emissivity will 
be somewhat larger (±7 K) due to the reduced zenith opacity. For unambiguous 
detection, the minimum detectable ice IWD at 90 GHz rises to ~ 0.02 g/cm 2 . Thus, 
even when water vapor and surface emissivity variations are considered, the SMMW 
sensitivity to ice IWD for polar winter conditions is approximately a factor of five 
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greater than at 90 GHz. 


For subtropical summer conditions (open circles), cloud water can be liquid, 
ice, or mixed phase, depending on altitude. For the 8-km base-altitude cloud, virtually 
all cloud water is frozen. However, for the 6-km base-altitude cloud, an appreciable 
fraction of the column (~ 55%) is liquid. At 90 GHz, the larger extinction coefficient 
of cloud liquid (relative to ice) causes the lower altitude (6 km, mixed-phase) cloud 
to be more easily detected than the high-altitude (8 km, ice-phase) one (Fig. 6a). 
Conversely, at higher frequencies (e.g., 166, 220, 340, and 410 GHz), two separate 
effects allow for easier detection of the higher-altitude cloud: (1) the single-scattering 
albedo of ice is significantly larger than liquid. Thus, the ice cloud top is more 
reflective than the ice-liquid cloud, and (2) the opacity of the increased water vapor 
burden between the satellite and the lower- altitude cloud top is greater. Thus, the 
response to the lower cloud is diminished. The latter effect is referred to as “water 
vapor screening”. Channels from 166 to 410 GHz are well suited to detecting high- 
altitude ice clouds due to an absence of water vapor screening. Lower-altitude liquid 
and mixed-phase clouds are also detectable at 166-410 GHz higher frequencies, but the 
advantage in sensitivity relative to 90 GHz is expected to be completely lost for cloud 
tops below ~ 4 -6 km from screening. The greatest sensitivity to subtropical summer 
clouds of arbitrary altitude is expected at ~ 220 GHz. The greatest invariance to 
cloud altitude (i.e., cloud water phase) is expected at ~ 166 GHz. 

Neglecting water vapor and surface emissivity variations, the minimum de- 
tectable IWD’s for high-altitude (8 km base altitude) subtropical-summer ice clouds 
are comparable to those found for the polar winter case. However', for unambiguous 
cloud detection in the presence of water vapor and surface variations, the subtropical- 
summer detection thresholds must be increased to ~ 10- 25 K, depending on the 
channel frequency. In this case, the minimum detectable IWD rises to ~ 0.1 - 
0.15 g/cm 2 for any of the channels from 90 to 410 GHz, and the sensitivity advan- 
tages of the high frequency channels (166-410 GHz) appear to be lost. However, 
these advantages can be recovered if the water vapor distribution can be indepen- 
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dently measured, as discussed in Sec. 6. Note that for subtropical summer conditions, 
the 340 and 410 GHz channels are highly insensitive to variations in either surface 

temperature or emissivity. 

The analysis over an ocean background modelled by a smooth Fresnel-refiecting 
surface yields similar results. Here, the surface reflectivity is 30-45% at 90-GHz, but 
decreases to 15-25% at 340-GHz. At the high EHF and SMMW frequencies, this is 
only slightly greater than for the case of land, hence the resulting cloud sensitivities 
are nearly identical. Any differences are even further reduced by ocean foam coverage. 
In contrast, the sensitivity to clouds at 90-GHz is significantly affected by the highly 
reflecting ocean surface. Not only is the brightness cooling caused by scattering ice 
clouds reduced, but a brightness warming caused by absorbing water clouds becomes 
observable. Although this bimodal effect complicates comparisons between oceanic 
cloud water sensitivities at 90 and 340 GHz, it appears that the sensitivity advantages 
of the SMMW channels that have been shown over land are only further enhanced 

over ocean. 

It is also recognized that random surface and atmospheric temperature vari- 
ations can produce brightness perturbations comparable to those caused by surface 
emissivity variations. However, the effects of surface and atmospheric temperature 
variations are expected to be more easily correctable than those of surface emissivity. 
Indeed, ocean surface temperature is a slowly varying process that can be determined 
to within ~ 1° C . by a combination of shipborne measurements, satellite measure- 
ments, and climatology. In addition, an increase in ocean temperature causes the 
Fresnel ocean emissivity at EHF and SMMW frequencies to decrease, thereby partly 
cancelling the effects of the temperature increase on the overall emission. Soil tem- 
peratures and vegetation-covered surfaces are highly variable on a diurnal time scale, 
but exhibit correlations with lower-tropospheric temperatures; both surface and tro- 
pospheric temperatures are observable to within ~ 1 - 3 K by passive microwave 
sensors operating at lower frequencies [24, 25]. For all of these surfaces, the impact 
of residual surface temperature variations on cloud detection capabilities should be 
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minimal. Variations in the temperatures of wet soil and vegetation-, snow- and ice- 
covered surfaces are slightly more difficult to characterize. However, upon correction, 
their impact on the cloud detection thresholds should not dominate over those of 
surface emissivity. Overall, the assumption of a maximum surface emissivity error of 
±3% justifies neglecting surface temperature variations when corrections are used. 

The preceding discussion assumes that the responses of the SMMW channels 
to IWD are independent of the actual thickness of the cloud. Indeed, for a given 
(fixed) hydrometeor size distribution, this assumption is well founded, as illustrated 
by the computed brightness perturbations for a subtropical summer cloud with a 
middle-altitude of 7.5 km, a constant IWD of 2.0 g/cm 2 , and thickness varying from 
0.5 to 4.0 km (Table 2, case “a”, using MP and SS mean sizes for 0.1 g/m 3 density). 
Here, the variation in the SMMW brightness perturbation remains below ±14% of its 
average as the thickness is varied, suggesting that IWD can be accurately measured. 
Of course, only the IWD within the top one to two optical depths can be directly 
probed. If the hydrometeor size distribution is allowed to vary with density (case “b”, 
using the MP and SS size relations), the resulting SMMW brightness perturbations 
vary over ±44% of their average. It is concluded that an independent measurement 
of mean hydrometeor particle size appears to be required for unambiguous IWD 
measurements. The potential for such a measurement using a widely spaced set of 
SMMW channels is discussed in Sec. 6. 

The cloud model used in this section assumes a linearly mixed phase between 
the freezing level and the altitude of ice nucleation (assumed to be -30 0 C ). While 
this is an accepted model for precipitating clouds, the mixed phase is unusual in 
non-precipitating clouds. Here, the Bergeron-Findeisen process rapidly depletes the 
liquid phase, resulting in precipitating ice that subsequently liquifies upon descent. 
However, insofar as passive sensing is concerned, the effects of a mixed-phase cloud 
are similar to those of a two-layer ice-over-liquid cloud. Such a two-layer cloud does 
not exhibit Bergeron formation, and is thus relatively stable. The reason for the 
radiometric similarity is that passive EHF and SMMW observations are not extremely 
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sensitive to the precise height of partially opaque layers of hydrometeors. Thus, for 
the purposes of this section, it is appropriate to use the linearly mixed-phase model. 
If the mixed-phase region were deleted from the model, then the proper choice of ice 
nucleation temperature becomes important. By choosing 0° C , the model becomes 
biased toward strongly-scattering ice clouds. By choosing -30° C , the model becomes 
biased toward absorbing supercooled liquid clouds. While uncommon in nature, the 
mixed-phase cloud model effectively sets the ice nucleation temperature in between 
these extremes. 


5 Precipitation Mapping 


Channels at 220, 340, and 410 GHz are more sensitive to clouds than the con- 
ventional EHF channels, although they are not as sensitive as IR channels, and are not 
expected to exhibit brightness temperature saturation for many cases of convective 
and stratiform precipitation. To demonstrate, brightness temperatures over simulated 
heavy clouds and raincells have been computed using the iterative numerical model 
for the high-frequency window channels 90, 166, 220, 340, and 410 GHz (Fig. 7). 
The computations used the iterative method, assuming MP liquid spheres, SS ice 
spheres, and HG phase functions with the total hydrometeor density uniform over 
the indicated altitude regions. The temperature profile is interpolated from the U.S. 
standard atmosphere listings for July at 35° N latitude, and the temperature of ice 
nucleation is -30° C . An exponential water vapor decay with a 2-km scale height and 
20% nominal surface relative humidity, and a land background with 5% reflectivity is 
assumed. The accuracy of the model has been verified for a case of mature convective 
precipitation using nadiral 118-GHz spectra [1]. The model also yields results that 
Me consistent with a numerical microwave radiative transfer model currently used at 
the NASA Goddard Space Flight Center [26]. 

Several scenarios are displayed in Fig. 7: moist and dry clear-air (curves 1 
and 2), and both low- and high-altitude clouds (curves 3, 4, 6, 7, and 9) and precip- 
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itation (curves 5,8, and 10). The hydrometeor density varies from 0.1 to 1.0g/m 3 , 
corresponding to MP precipitation rates of ~ 1 - 18 mm/hr, respectively. As shown 
by Adler et al. [27] using coincident microwave and IR data, an IR channel would be 
expected to saturate over any of the hydrometeor-laden scenarios (3-10) in Fig. 7 (see 
also [28]). Conversely, microwave channels at 92 and 183 GHz have been shown 
to contain useful information on the underlying storm structure and ram rate (see, 

e.g., [29]). 

From Fig. 7, the response of the SMMW channels (340 and 410 GHz) relative 
to the conventional channels can be predicted. Over the eight cloud and raincell 
scenarios, the expected brightness variation of a SMMW channel is approximately 
one-third to three-fourths that of a 90-GHz channel. Further, the sensitivity of a 
SMMW channel to precipitation will be suppressed relative to 90-GHz, and the sen- 
sitivity to thin (non-precipitating) clouds will be enhanced. Overall, the SMMW 
channels will be somewhat less sensitive to storm structure and rain rate than at 
90-GHz, but more sensitive to the presence of ice canopies. 

The computed weighting functions for absorbing and scattering hydrometeor- 
laden atmospheres at 90- and 340-GHz (Fig. 8a and b) help explain these sensitivity 
differences. Two features are essential in determining the upwelling brightness: the 
weighting function peak-altitude and the cosmic background weight W C b (the back- 
ground weight is essentially equal to the cell top reflectivity since little reflection 
occurs via the land surface). At low ice densities (< 0.3 g/m 3 , or ~ 1 mm/hr SS ice 
precipitation rate), scattering and absorption at 340 GHz (Fig. 2b) cause up to 20% 
cell top reflection, along with a high-altitude (~ 8 — 10 km) weighting function peak. 
Both of these effects produce a cooling in the upwelling brightness, although the mech- 
anism of cell top reflection, rather than emission from colder levels, is dominant. At 
higher ice densities, the 340-GHz single-scattering albedo decreases slightly, effectively 
saturating the cell top reflectivity. Any further decreases in 340-GHz brightness are 
due primarily to increases in the cell-top opacity, which further attenuate radiation 
originating from warm low altitudes. 
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At 90 GHz, ice scattering is insignificant until a critical density of ~ 0.5g/m 
is exceeded. Such a critical value arises due to a strong nonlinearity between cloud 
density and the volume scattering coefficient: as cloud density increases, so does the 
mean cloud particle size. For Rayleigh scatterers, this causes an abrupt increase 
in the scattering coefficient (as illustrated by Fig. 5b). Thus, the 90-GHz brightness 
exhibits little change at lower densities. At higher densities, the cell is radiometrically 
opaque [1]: the weighting function peaks within the ice layer, and the reflectivity 
rapidly increases with density. However, in this case the 90-GHz ice single- scattering 
albedo does not significantly decrease with increasing density. Thus, at densities 
greater than ~ 0.5 g/m 3 , the 90-GHz cell top reflectivity exceeds that at 340 GHz 
by several percent, and the resulting brightness temperature does not saturate for 
typical cell densities. 

Screening by water vapor above the cell top contributes further to reduced 
sensitivity at 340 GHz relative to 90 GHz. By considering the summer low-latitude 
atmospheric profile in Fig. 7, the extreme effects of water vapor screening on SMMVV 
observations are displayed. For example, curves 3-5 show the window-channel re- 
sponse to low-altitude (6 km cell top) liquid precipitation. Although the liquid 
drops exhibit significant scattering and absorption at SMMW frequencies, variations 
in brightness are screened by intervening water vapor. For high-altitude cells, the 
screening opacity is small (curves 6-8,10), and the resulting variations in the SMMW 
channels are significantly greater. For these scenarios, the SMMW channels will not 
saturate, and can provide useful information on storm structure, heavy clouds, and 
precipitation for many scenarios in which IR channels would saturate. Precipitation 
cell cores should be detectable at 220-, 340-, and 410-GHz beneath many optically 
opaque (0.01-0.05 g/cm 2 ) cirrus canopies. By virtue of their potential for enhanced 
spatial resolution, such channels would be particularly useful for precipitation map- 
ping and resolving retrieval ambiguities caused by fractional footprint-filling in 10-, 
18-, and 37-GHz images of oceanic rain rate. At high-latitudes or during winter con- 
ditions, the screening opacity due to water vapor will be considerably lower, and the 
SMMW channels will be particularly useful for detecting and mapping both high- 
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and low-altitude precipitating clouds. 


For radiometrically opaque cells (i.e., Ws « 0), the cell probing depth is the 
distance down into the cell over which most emission occurs. This defines a region 
which comprises most of the weighting function’s area below the cell top. Below this 
region, no direct radiometric observations of hydrometeor parameters can be made. 
Under most conditions, the SMMW channels will not be able to directly probe sur- 
face rain rate. However, indirect information on low-level rainrate might be inferred 
statistically from correlations with direct radiometric observations of other ramcell 
parameters, e.g., the cell’s size and altitude [30, 31]. In addition, the monotonic 
decrease of cell probing depth with increasing frequency suggests that hydrometeor 
parameter profiling in heavy clouds and precipitation cell tops might be facilitated 
using conventional channels along with a widely spaced set of SMMW channels. The 
concept of passive precipitation profiling has been discussed by Kummerow [29]. 

6 Additional Meteorological Degrees of Freedom 
through Wideband Observations 


The enhanced sensitivity to ice clouds at 220 and 325 GHz is a deterrent 
for water vapor profiling. However, it is expected that by using these channels in 
conjunction with the conventional channels (90, 166, and 183 GHz), additional me- 
teorological information (not obtainable using either of these bands exclusively) can 

be obtained. 

Support for this notion arises from nadiral aircraft observations at 60 and 
118 GHz [1, 32]. The data has revealed that perturbations caused by clouds and light 
precipitation are typically 2-3 times as large in the 118-GHz channels as compared to 
corresponding 60-GHz channels (i.e., those with similar clear-air weighting functions). 
The observed sensitivity differences have been corroborated by calculations using the 
iterative radiative transfer model, and are due to Rayleigh scattering and absorption 
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by the electrically small hydrometeors. Thus, coincident observations of clouds and 
light precipitation using similar 60 and 118-GHz channels results in an additional 
degree of freedom related to liquid and ice water content, although it is not clear at 
this time how to best use this information. Over heavy precipitation, the scattering 
and absorbing hydrometeors are electrically large. In this case the coincident 60 and 
118-GHz perturbations are comparable, and provide no additional information. 

Analogously, at least one additional degree of freedom is expected using sim- 
ilar channels near the 183- and 325-GHz HjO lines. Due to the Rayleigh frequency 
dependence of small (< 0.1 mm diameter) hydrometeors, these channels will respond 
differently to layers of non-precipitating cloud. Referring again to Fig. 4, a thin ice 
cloud layer (8-10 km, 0.1 g/m 3 ) perturbs channels near 325 GHz more strongly than 
near 183 GHz (nearly a 15 K difference), suggesting that information concerning the 
presence of the thin cloud is available from the difference spectra between these two 
bands. This information can potentially be used to correct water vapor soundings, 
and is only available from coincident observations at both 183 and 325 GHz. As the 
cloud density increases (8-10 km, 0.3 g/m 3 cloud layer and 1-10 km, 1 g/m 3 cell), 
the sign of the trend in Fig. 4 reverses, and the 183-GHz channels are more strongly 
perturbed (nearly a -20 K difference). The bimodal sensitivity is a consequence of the 
transition from Rayleigh to Mie scattering as the hydrometeor density (and hence size) 
increases. For frequencies below the Rayleigh-Mie transition peak (Figs. 2a and b), 
scattering and absorption increase with frequency. Above this peak, scattering by ice 
and both scattering and absorption by liquid decrease slowly with frequency. 

Another example of an additional degree of freedom is illustrated in Fig 9. 
Computed brightness perturbation spectra over a rain cell using a wideband channel 
set are shown for two different ice size distributions: MP and SS. The model cell 
density is 1 g/m 3 and extends from 1 to 10 km in altitude. The temperature of ice 
nucleation is -30° C. The MP-distributed ice spheres (which are typically smaller than 
SS-distributed spheres for densities near 1 g/m 3 ) cool the brightness spectrum at the 
higher EHF frequencies (90, 118, 166, 183, 220 GHz) and increase the brightness 
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spectrum at lower frequencies (6, 10, 18, 22, 37 and 60 GHz). Thus, the brightness 
perturbation spectrum “pivots” around a frequency near 75 GHz. The slope of this 
spectrum, as measured using a widely spaced channel set, contains information on 
the particle size distribution. 

For a wideband EHF and SMMW channel set, a similar degree of freedom 
related to cloud ice particle size might be expected. This is illustrated in Fig. 7 
(curves 7 and 9), for a 2-km thick, 0.3 g/m 3 ice layer. Here the computed window 
channel spectrum also pivots as the ice size distribution is changed from SS to MP, 
but around a frequency of ~ 140 GHz. The change reduces the mean ice radius 
from ~ 0.22 mm to ~ 0.16 mm, which subsequently increases the 166, 220, 340, 
and 410 GHz scattering coefficients, thereby cooling the high-frequency portion of 
the spectrum. Simultaneously, scattering is slightly decreased at 90 GHz, thereby 
warming the low-frequency portion. 


Generalizing, we can define the spectral pivot frequency f p as the frequency at 
which the hydrometeor extinction exhibits minimum sensitivity to mean hydrometeor 
size: 

: 0 (5) 




d(a) 

From the preceding discussion, it appears that f p fts 32/(a), where a is in millimeters 
and f p is in gigahertz. This is in general agreement with relationships determined 
graphically using constant-extinction contours (Figs. 10a and b), where it is seen 
(dashed lines) that f p « 15/(a) for liquid and f p « 25 /(a) for ice. In either case, f p is 
slightly below the Rayleigh-Mie transition frequency, and results from the scattering 
peak near this frequency. 


The sensitivity to size distribution displayed by the perturbation spectra in 
Figs. 9 and 7 suggest an observable mode which is related to the mean ice particle 
size. Indeed, the slope of the brightness spectrum near f p is sensitive to the mean 
hydrometeor size, and thus might potentially be useful for such measurements. Due 
to the inherent lack of structure in microwave hydrometeor scattering and absorp- 
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tion spectra, a widely-spaced window-channel Bet (preferably spanning at least one 
octave) with high spatial resolution (to eliminate footprint-filling ambiguities) would 
be required to observe this hypothesized hydrometeor size degree of freedom. The 
EHF and SMMW window channels from 90 to 410 GHz are well suited for this task. 

Brightnesses near f v exhibit minimum sensitivity to hydrometeor size varia- 
tions, and thus should be useful for measurements of other cell parameters such as 
IWD. For many non-precipitating clouds, f p is between ~ 200 — 600 GHz, suggest- 
ing that the SMMW channels 220, 340 and 410 GHz are particularly well suited for 
non-precipitating cloud IWD measurements. 


7 The Millimeter- Wave Imaging Radiometer (MIR) 


The Millimeter- Wave Imaging Radiometer (MIR) is a new instrument be- 
ing designed for studies of airborne passive microwave retrieval of tropospheric wa- 
ter vapor, clouds, and precipitation parameters. The MIR is a total-power cross- 
track scanning radiometer for use on either the NASA ER-2 (high-altitude) or DC-8 
(medium altitude) aircraft. The current design includes EHF channels at 90, 150, 
183±1,3,7, and 220 GHz, and SMMW channels at 325±1,3, and 9 GHz. Identical 
3.5° beamwidths for all channels produce nadiral 3-dB footprint sizes of ~ 700 m at 
mid-altitude (~ 10 km) when operated aboard the NASA ER-2 (the cruise altitude 
of the ER-2 is ~ 20 km). The ±50° field-of-view is scanned every ~ 2.5 seconds. 

The current MIR design consists of a scanhead and data acquisition system, 
designed for installation in the ER-2 superpod nose cone. The scanhead will house 
the receivers (feedhorns, mixers, local oscillators, IF amplifiers, and video detectors), 
a scanning mirror, hot and cold calibration loads, and temperature sensors. The 
data acquisition system is to consist of a 80286-based computer with 1000 MB tape 
drive, a 12-bit 48-channel opto-isolated A/D converter, 8-bit D/A converter (for offset 
drift compensation) and a power supply /conditioning module. Software-controlled 
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offset and gain adjustments will compensate for expected receiver gain and noise 
temperature fluctuations. A watch-dog timer and IRIG-B time decoder are included. 

The high spatial resolution EHF and SMMW brightness imagery will be useful 
for: (1) Verifying EHF and SMMW radiative transfer models in clear air, clouds and 
precipitation, with direct comparisons at frequencies up to 334 GHz, (2) Constructing 
and evaluating EHF and SMMW mapping and retrieval algorithms for precipitation 
parameters (e.g., rainfall rate, particle size, and cell-top altitude), cloud parameters 
(e.g., liquid and ice water content, and particle 6ize), and water vapor profile retrievals. 
Of particular interest are potential measurements of cirrus ice water content, and eval- 
uation of the 325-GHz water vapor line for spaceborne profiling of tropospheric water 
vapor. In addition, the imagery will be useful for: (3) Studying convective raincell 
and atmospheric gravitational wave structure, (4) Evaluation of passive microwave 
satellite data from instruments such as AMSU-B and SSM/T2, and (5) Develop- 
ment of techniques for combining observations at EHF and SMMW channels along 
with IR and conventional microwave channels to facilitate the estimation of temper- 
ature and water vapor profiles, and precipitation parameters, particularly ice density 
and ice size distributions. Many of the MIR science objectives will require com- 
plementary data from other airborne, ground-, and satellite-based instruments (e.g., 
imaging radiometers operating at complementary frequencies, ground-based and air- 
borne doppler weather radars, satellite- and airborne optical and infrared imagers, 
and radiosondes). 


8 Discussion 


In this paper, we have outlined the chief advantages and disadvantages of pas- 
sive EHF and SMMW channels as applied to satellite-based tropospheric and lower 
stratospheric remote sensing. Briefly, the following conclusions are drawn concerning 
these largely unexplored channels: (1) Enhanced spatial resolution from diffraction- 
limited antennas of practical size is available using the SMMW channels. (2) Channels 
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around the 325-GHz water vapor line and in the 220-GHz window can be used for 
water vapor profiling in clear air with an accuracy comparable to that expected using 
166 and 183-GHz channels. However, with no additional information, the adverse 
effects of clouds will be greater. (3) Clouds can be detected using 340 or 410 GHz at 
densities as low as 0.004g/cm 2 , particularly under polar winter conditions or at alti- 
tudes greater than ~ 6 km. The available sensitivity to clouds increases significantly 
with frequency above 90 GHz. Thus, the SMMW channels are expected to be useful 
for imaging cirrus and non-precipitating cumulus structure. For subtropical summer 
clouds under ~ 4 - 6 km altitude, 90, 166, and 220 GHz appear best for detection 
and IWD measurement. (4) Many heavy clouds and precipitation can be mapped 
using SMMW channels with one-third to three-fourths of the brightness variation ex- 
pected using 90 GHz. Saturation, as exhibited in IR imagery of rain cells, is not often 
expected. Images of cirrus clouds using SMMW channels and raincells might poten- 
tially facilitate the indirect retrieval of surface rain rate using other low-frequency 
microwave channels (e.g., 10, 18 and 37 GHz). However, contrast reduction due to 
water vapor screening will be significant for cells with low top-altitudes (~ 6 km or 
less) during subtropical summers. (5) Additional meteorological degrees of freedom 
related to water vapor, cloud integrated water density, and mean particle size are 
expected using a wideband set of EHF and SMMW channels. For unambiguous de- 
tection of clouds under summer subtropical conditions, independent information on 
the water vapor profile will be required. This could potentially be obtained from such 
a wideband set. It also appears possible that a wideband set of EHF and SMMW 
channels (e.g., 90, 166, 220, 340, and 410 GHz) can be used to retrieve both cloud 
integrated water density and mean hydrometeor size. 

Major deterrents to developing meteorological retrieval algorithms using EHF 
and SMMW passive channels are the nascent states of radiative transfer science and 
the lack of experimental radiometric imagery for microwave frequencies above 200 
GHz. Meteorological data flights of the MIR aboard the NASA ER-2 along with sub- 
sequent data analyses will be crucial in refining radiative transfer models (particularly 
for cirrus clouds) and in developing EHF and SMMW retrieval algorithms. The first 
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MIR data flights are tentatively planned for 1992. Measurements of mean ice par- 
ticle size and cloud IWD could provide useful information on the infrared radiative 
properties of cirrus clouds; these results could have an impact on the role of cirrus in 
global warming and links between cirrus clouds and precipitation. Currently, cirrus 
cloud radiative properties are a major source of uncertainty in global radiation and 

climate models. 

We specifically have not addressed the issues of nascent SMMW receiver and 
antenna technology, or the increased integration noise expected for high-resolution 
SMMW observations due to reduced spot dwell time. Indeed, enhanced spatial res- 
olution using a single-beam scanning imager can be gained only at the expense of 
reduced spot integration time, and hence increased radiometric noise. However, such 
increases can potentially be offset by using multiple scanning beams, wider IF band- 
widths, longer scene revisit times, and improvements in receiver sensitivity. Reflector 
surface accuracies of 10/im or less are acceptable up to 410 GHz, and are technolog- 
ically feasible for spaceborne antennas several meters in diameter. It is (somewhat 
tacitly) assumed that these and other associated problems are surmountable through 
future improvements in radiometer and large space structure technology. However, 
careful absolute calibration and high radiometric resolution will be required to observe 
any wideband spectral modes. 
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Table 1: Line parameters b lh bn, hi and center frequencies 1 /, for seven of the most sig- 
nificant H 2 0 lines for tropospheric remote sensing, along with continuum parameters 


bj and b e [12]. 


Vi 

(GHz) 

bu 

(kHz/mbar) 

b 2i 

hi 

(Ghz/mbar) 

22.2351 

0.0109 

2.143 

0.002784 

183.3101 

0.23 

0.653 

0.003164 

325.1529 

0.154 

1.515 

0.00297 

380.1974 

1.19 

1.018 

0.003036 

448.0011 

1.06 

1.370 

0.00238 

556.9360 

51 

0.114 

0.0030 

752.0332 

25 

0.336 

0.00286 

bj = 

1.13 x 10' b 

(GHz-’mbar-*) 

h = 

3.57 x 10" 7 

(GHz -1 mbar -2 ) 


Table 2: Brightness perturbations for a 7.5-km altitude subtropical summer cloud of 
varying thickness and fixed integrated water density (liquid and ice) of 2-0 g/cm . 
Calculations assume mean hydrometeor sizes that are (a) fixed, at the MP/SS values 
for 0.1 g/m 3 density, and (b) varying with density, according to the MP/SS relations. 


Size 

Distribution 

Mean Liquid 
Radius 
(mm) 

Mean Ice 
Radius 
(mm) 

Cloud 

Thickness 

(km) 

Brightness Decrease (K) 
at Frequency (GHz): 

90 166 220 340 410 


Clear-air: no 

cloud 


0 

0 0 0 u 

(a) 

( a ) 

(a) 

(a) 

(a) 

0.125 

0.126 

0.5 

-14.3 

-26.7 -30.5 -25.3 -18.5 

0.125 

0.126 

1.0 

-14.5 

-27.2 -31.2 -25.9 -19.1 

0.125 

0.126 

2.0 

-14.7 

-28.2 -32.4 -27.4 -20.8 

0.125 

0.126 

3.0 

-15.0 

-29.2 -33.7 -29.2 -22.6 

0.125 

0.126 

4.0 

-15.1 

-30.1 -34.8 -30.9 -24.4 

(b) 

(b) 

(b) 

(b) 

(b) 

0.177 

0.261 

0.5 

-16.1 

-21.8 -21.1 -14.4 -10.2 

0.149 

0.182 

1.0 

-16.0 

-25.9 -27.2 -20.2 -14.7 

0.125 

0.126 

2.0 

-14.7 

-28.2 -32.4 -27.4 -20.8 

0.113 

0.102 

3.0 

-13.7 

-28.7 -34.5 -31.9 -25.1 

0.105 

. . « « 1 • • 

0.088 
n J _ a % 4 

4.0 

rn /cc 

-12.8 

. n 1 

-28.6 -35.5 -35.1 -28.7 


(a) Mean liquid and ice radii fixed at MP/SS values for 0.1g/m 3 density. 

(b) MP-distributed liquid and SS-distributed ice spheres. 
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Figure 1: Tropospheric clear-air EHF and SMMW absorption for several relative 
humidities. Specific frequencies of interest are indicated by arrows. 


Figure 2: Polydispersive Mie hydrometeor absorption and scattering (from Gasiewski 
and Staelin [1]): (a) Liquid, assuming a Marshall-Palmer drop size distribution. Com- 
putations from Savage [2] for absorption (•) and scattering (x) are plotted for com- 
parison. (b) Ice, assuming a Sekhon-Srivastava distribution. Calculations are shown 
for precipitation rates of 1, 10, and 40 mm/hr for both phases and 100 mm/hr for 
liquid. 


Figure 3: Clear-air temperature weighting functions vs. altitude for channels at 
(a) 90, 166, and 183±1,3,7 GHz, and (b) 220, and 325±1,3,9 GHz. The right-hand 
ordinate is the approximate water vapor burden. 


Figure 4: Computed nadiral brightness spectra differences between similar 183- and 
325-GHz channels for clear-air (20% and 80% surface relative humidity), clouds (0.1 
and 0.3 g/m 3 density, 8-10 km altitude), and precipitation (1 g/m 3 density, 1-10 km 
altitude). 


Figure 5: Polydispersive Mie hydrometeor absorption and scattering for a 0.1 g/m 3 
cloud of (a) liquid and (b) ice, assuming exponential size distributions. Calculations 
are shown for average particle radii (rad) of 0.01, 0.03, 0.1 and 0.3 mm at 0 C . 


Figure 6: Computed nadiral brightness changes due to non-precipitating clouds 
(solid lines) and water vapor (dashed lines) at (a) 90, (b) 166, (c) 220, (d) 340 and 
(e) 410 GHz. Open circles are for a subtropical summer atmosphere (July, 30° N 
latitude) and dark circles are for a polar winter atmosphere (January, 90° N lati- 
tude). Clouds are 2 km thick with base altitudes at either 6 or 8 km. The abscissa is 
either cloud density or surface relative humidity (as applicable). (Note the different 
ordinate scales.) 
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Figure 7: Computed window-channels nadiral brightnesses over heavy clouds and 
convective precipitation: 

(1) clear-air, 20% surface relative humidity, 

(2) clear-air, 80% surface relative humidity, 

(3) 4-6 km, 0.1 g/m 3 , MP liquid cloud layer, 

(4) 4-6 km, 0.3 g/m 3 , MP liquid cloud layer, 

(5) 1-6 km, 1 g/m 3 MP-SS cell, 

(6) 8-10 km, 0.1 g/m 3 , SS ice cloud layer, 

(7) 8-10 km, 0.3 g/m 3 , SS ice cloud layer, 

(8) 1-10 km, 1 g/m 3 MP-SS cell, 

(9) 8-10 km, 0.3 g/m 3 , MP ice cloud layer, 

(10) 1-14 km, 1 g/m 3 MP-SS cell. 


Figure 8: Computed nadiral temperature weighting functions at (a) 90 GHz and 
(b) 340 GHz for clear-air, clouds, and precipitation. The symbols near 0 and 25 km 
indicate the surface and cosmic background contributions (respectively). 


Figure 9: Nadiral brightness perturbation spectra for a model 1-10 km, 1.0 g/m 3 
precipitation cell. The calculations were made using the iterative planar- stratified 
numerical model for two ice particle size distributions: Marshall-Palmer (MP) and 
Sekhon-Srivastava (SS). 


Figure 10: Extinction contours for exponentially-distributed Mie (a) liquid and (b) ice 
spheres. Contour values are in dB-km“ 1 -g“ 1 -m 3 . Along the dashed lines, the sensi- 
tivity of extinction to variations in mean hydrometeor size vanishes. 
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Appendix B. 


(Proceedings of the 1992 tree Joint Symposia APS/URSI/IMP 
Meeting, Chicago, Illinois, July 18-25, 1992.) 

Electromagnetic Scattering from Microwave 

Absorbers: 

Laboratory Verification of the Coupled Wave Theory 


A.J. Gasiewski, D.M. Jackson* 
School of Electrical Engineering 
Georgia Institute of Technology 
Atlanta, GA 30332-0250 


Laboratory W-band measurements of the bistatic scattering func- 
tion of some common microwave absorbing structures have been made. The 
structures investigated include periodic wedge-type and pyramid-type iron- 
epoxy calibration loads and flat carbon-foam “Echosorb” samples (Table 1). 
Measurements were made using an HP 8510 network analyzer interfaced to 
a focussed-lens scattering range (Figure 1). The test set provided a dynamic 
range of ~ 55 dB. Swept frequency measurements over the band 75-100 GHz 
have revealed specular and Bragg reflection characteristics in the measured 

data. 

The range consists of two mechanical booms (one each for the re- 
ceiver and transmitter) that rotate about a common axis. A standard-gain 
horn and rexolite lens (focal length = 8.5 cm) are mounted on each boom 
and are aligned to produce a 4-cm (approximately) focussed spot on the sam- 
ple under test, located on the axis of rotation. Both co- and cross-polarized 
reflection measurements for both TE- and TM-incident polarizations can be 
made. Figure 2 (A1 plate) shows the 90-GHz angular response of the range 
for a fixed incident angle (60°) and variable scattering angle. 

In interpreting the data, the angular response of the system to the 
specular scattering from a flat aluminum plate is deconvolved from the mea- 
sured reflectivities. Three deconvolution procedures are considered: Fourier 
transformation, linear minimum mean-square error estimation (LMMSE), 
and the CLEAN algorithm [1]. A major drawback of Fourier transforma- 
tion is an enhancement of noise energy near the zeroes of the transformed 
angular response. The LMMSE technique is an implementation of a dis- 
crete non-causal Wiener filter and is less sensitive to measurement noise. 
The CLEAN algorithm was originally developed for radioastronomical im- 
age enhancement in which the measured data consisted of a two-dimensional 
distribution of point sources convolved with the response function of the ra- 
diotelescope. Since the bistatic scattering function of a periodic surface is an 
angular distribution of impulses, CLEAN is also well-suited for deconvolution 
of the W-band scattering measurements. 

Theoretical calculations of the bistatic scattering coefficient for the 
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periodic wedge samples were performed using the coupled wave method [2, 3]. 
In this method, the periodic region between the peaks and troughs of the 
wedge structure is discretized into inhomogeneous horizontal layers. The 
permeability and permittivity in the layers are each represented as Fourier 
series. The EM fields within each layer axe expressed as a sum of inho- 
mogeneous plane waves. The fields above and below the wedge region are 
expressed as a sum of propagating and evanescent Floquet harmonics. By 
using Maxwell’s curl equations, a set of coupled equations involving the un- 
known wave amplitudes is obtained. A solution for the amplitudes is foun 
using eigenanalysis and EM-field boundary conditions at the interfaces to 
each layer. The numerical accuracy for the fields in each region is checked 
using conservation of energy. 

Figure 2 is a comparison between calculated and measured co- 
polarized TM reflection at 90-GHz for sample #4 (wedge absorber), assum- 
ing tr = 9.0 - j'0.4 and fi r = 1.0 - j0.5. The computed location of the 
first Bragg lobe (-10.1°) is consistent with the measured response near this 
angle. Moreover, the calculated magnitudes of both the specular lobe and 
first Bragg lobe are consistent with the measured values. For this case, all 
other reflected harmonics are evanescent. The measured reflectivity for sam- 
ple #5 showed evidence of several propagating harmonics, as predicted. For 
example, at 90 Ghz and incident angle of 60°, there are 7 Bragg lobes. The 
laboratory measurements for the wedge absorbers also showed higher specu- 
lar and Bragg reflection for TE-incident waves than for TM-incident waves, 
as expected. 

Additional measurements were made for several other structures. 
Measured TM reactivities for sample #2 (flat dielectric) showed only a small 
Brewster-angle null (~3 dB), indicating a complex permeability. Sample #3 
(flat open-cell absorber sheet) exhibited a TE co-polarized specular reflectiv- 
ity of ~ -10 dB at 90-GHz for 45° incidence. Due to the large electrical size 
of the cells, diffuse scattering was expected, but not observed in the raw data. 
For sample #6 (pyramid load), the dynamic range was insufficient to char- 
acterize the reflectivity. This was not unexpected since two-dimensionally 
periodic structures generally scatter energy over the entire upper hemisphere. 
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Sample # 

Description 

1 

Flat A1 plate 

2 

Flat Fe-epoxy slab 

3 

Flat carbon -foam “Echosorb” sheet 
(pore size ~ 0.5 mm) 

4 

Periodic wedge Fe-epoxy absorber 
(pitch=3.2 mm, depth=6.5 mm) 

5 

Periodic wedge Fe-epoxy absorber 
(pitch=12.5 mm, depth=24.0 mm) 

6 

Periodic pyramid Fe-epoxy absorber 
(pitch=10 mm, depth=40 mm) 


Table 1: Desci 
measurements . 
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Figure 1: Schematic diagram of the focussed scattering range. 



Figure 2: Measured reflectivity at 90 GHz and 0 t = 60° for Sample #1 (A1 
plate) and #4 (wedge load) for TM-incident polarization. The grating vector 
is in the plane of incidence. The coupled-wave calculations for the specular 
and first Bragg directions (scattering angles 60° and -10.1°, respectively) 
are shown (A) for comparison. 
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